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Abstract 

We investigate time evolution of prepared vibrational state (system) coupled to a reservoir with 
dense spectrum of its vibrational states. We assume that the reservoir has an equidistant spectrum, 
and the system - reservoir coupling matrix elements are independent of the reservoir states. The 
analytical solution manifests three regimes of the evolution for the system: (I) weakly damped 
oscillations; (II) multicomponent Loschmidt echo in recurrence cycles; (III) overlapping recurrence 
cycles. We find the characteristic critical values of the system - reservoir coupling constant for 
the transitions between these regimes. Stochastic dynamics occurs in the regime (III) due to 
inevoidably in any real system coarse graining of time or energy measurements, or initial condition 
uncertainty. At any finite accuracy one can always find the cycle number k c when dynamics of 
the system for k > k c can not be determined uniquely from the spectrum, and in this sense long 
time system evolution becomes chaotic. Even though a specific toy model is investigated here, 
when properly interpreted it yields quite reasonable description for a variety of physically relevant 
phenomena, such as complex vibrational dynamics of nano-particles, with characteristic inter-level 
spacing of the order of 10 cm" 1 , observed by sub-picosecond spectroscopy methods. 

PACS numbers: 05.45.-a, 72.10.-d 
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Experimental studies of vibrational dynamics in various systems (ranging from relatively 
small molecules in a gas phase (near dissociation boundary) or water molecules clusters con- 
fined to interface, through nano-particles or large photochromic molecules) reveal in a time 
domain 10~ 13 — 10~ n s very rich variety of regimes Q], 0|. Observed in such systems 
seemingly irregular damped oscillation regimes can not be explained theoretically in the 
frame work of widely used models with reservoirs possessing continuous spectra j^], Q], 
Indeed in the case of a system coupled to the continuous spectrum reservoir, only smooth 
crossover between coherent oscillations and exponential decay is possible upon increasing of 
the coupling. Moreover generic complex dynamics is observed in the systems with charac- 
teristic inter-level spacing of the order of 10cm _1 , when the recurrence cycle period is in the 
range of sub-picoseconds, where the measurements are performed. The dense but 

discrete spectra with characteristic inter-level spacings of this order are typical for a wide 
variety of nano-particles with 10 2 — 10 3 degrees of freedom. In such cases one should expect 
recurrence cycles and much more rich and sophisticated time evolution. Motivated by these 
observations, our intent here is to examine joint system- reservoir evolution, i.e., recurrency 
cycles, when the energy is flowing back from the reservoir to the system. 

We investigate a simple (but yet non-trivial) model of a single level system coupled to a 
reservoir with discrete spectrum 

H = e° s btb s + £ e°XK + E Cn(bjb n + 6,6+) , (1) 

n n 

where e°, and e° are bare (i.e., without interaction) eigen states of the system and of the 
reservoir 6 + are corresponding creation operators), and C n are interaction matrix ele- 
ments. The Hamiltonian (Op) has non-zero diagonal elements and only one non-zero row and 
column. The secular equation to find the eigen values e of the Hamiltonian reads as 

*X<0 = e-E7% = O. (2) 

fc fc n 

where the energy is measured in units of mean reservoir inter-level spacing, and counted 
from the unperturbed system energy e° s . 

Time dependent amplitudes satisfy to the corresponding Heisenberg equations of motion 

i&s = E C n a n i if^n = C n a s + e^a n , (3) 

n 

supplemented by the initial condition 

a s (0) = 1; a„(0) = 0. (4) 



The solution of these equations of motion can be represented as 

x . expjiet) 

as{t) = ^llFW 1 ^ (5) 
where Eq. (jSJ) is the sum over the residues in the simple poles e = e n which are the roots of 
the secular equation (|2"1) . Because of level repulsion phenomenon, the dense spectra can be 
grouped into series of approximately equidistant levels [8j . If the system- reservoir dynamics 
is dominated by a single among these series, the Hamiltonian ([1]) can be simplified following 
an exactly solvable model proposed long ago by R.Zwanzig [9( of a system coupled (inde- 
pendent of the reservoir states, i.e. C n = C) to a reservoir with equidistant spectrum (i.e., 
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in our notation e° = n) . In fact the Zwanzig model treats a simplified version of well 
Caldeira-Legget Hamiltonian widely used in condensed matter physics and chemistry 
What is lacking, as far as we know, is an investigation of the recurrence cycle dynamics for 
the Caldeira-Legget Hamiltonian. Indeed in the standard approach the reservoir spectrum 
is assumed to be so dense that any recurrence is irrelevant on a characteristic measurement 
time scale. However it is not always the case and many systems (see [l| - {4]) lie in the 
intermediate range of the parameters, with discrete but dense spectrum of final states. For 
the ease of notation let us assume also that the level n = of the reservoir is in the resonance 
with the system energy level e s . Thus we get from Eqs. ([2]) and (jSJ) 

e = nC 2 cot(ne) , (6) 

and 

Mt)=2 £i+*»™+l/c)»- (7) 

Utilizing the Poisson summation formula one can replace the series terms in Eq. (J7J) by 
their Laplace transforms, as 



a s(t) = / 7^ — 2^2 i t — Tnvi wpi-^ikn + ie n t) (8) 



-00 1 + tt 2 C 2 + (e n /C) 2 
£ loo ^2 T (e/C) 2 -M^{e - 0) + iet] 



tt 2 C 2 + (e/C) 

where we replace the integration variable n — > e n and introduce the phase = e n — n which 
can be taken in the interval (0, 1/2). According to the eigen- value equation ([6]), we have 

de n n 2 C 2 + (e n /C) 2 



dn 1 + n 2 C 2 + (e n /C) 2 



(9) 



and 

e ^- 2 ^ = tn^- (10) 

Since t > the sum in the (JED includes only a finite number of the cycle amplitudes < k < 
[t/(27r)] (where [x] stands for the integer part of x). The integrand has two poles e = iwC 2 , 
and only the positive pole contributes into the integral. Furthermore, the k < cycles do 
not contribute into a s (t), and also for k > [t/2ir] the pole inC 2 lies outside the integration 
contour in the upper half plane. Combining everything we end up with 

[t/2-rr] 

where 

c r ^t^P "''-"-'- (12) 
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The integral 
(where > 1) 



2J) can be expressed in terms of the generalized Laguerre polynomials L\ 



ai k \t) = 2^4_i(2r fc ) exp (-r fc ) 9{r k ) , (13) 



where 

Tfc = n C 2 (t - 2nk) (14) 

is the local time for the fc-th cycle. The step function 8{x) enters Eq. (Fl3|) because a[ k \t) = 0, 
for t < 2A;7r. Note that the Poisson summation formula replaces the discrete series of the 
poles along the real axis by a single pole on the imaginary axis. This pole determines decay 

n 

probability of the initial quasistationary state [111 ]. The expressions ( TTTj) - f[T3"j) describe the 
time evolution of the system. The equation (ITTj) represents dynamics in terms of a certain 
superposition of the coherent eigen-frequency oscillations, whereas the equation (TlBl - in 
terms of the transition probabilities. 

In the limit of weak coupling, nC 2 1, the system time evolution is dominated by the 
coherent transitions between the resonance states (i.e., n = in our case). There are many 
recurrence cycles within period 2n/C of these oscillations. For the strong system - reservoir 
coupling, ttC 2 ^> 1 the transitions to many reservoir states, |n| < ttC 2 , contributes to the 
system time evolution. Interference between the transitions suppresses probability back-flow 
(from the reservoir to the system). In the initial cycle fc = 0we get exponential decay 

af(t) = e- r ', (15) 



where as one could expect, the exponential decay rate is determined by the Fermi golden 
rule 

r = vrC 2 . (16) 

However (and it is one of our new observations in this paper) in the following recurrence 
cycles we find a sort of Loschmidt echo. This phenomenom occurs due to quantum mechan- 
ical synchronization of the reservoir - system transitions. The Loschmidt echo, we found, 
occurs because in the strong coupling limit, back and forth transitions of many reservoir 
states which determine system dynamics, take place not at the same time. We illustrate 
these fine structure features of the Loschmidt echo on the Fig. 1, where a s (t) is calculated 
from Eq. f[T3"j) for C 2 = 0.3 and for C 2 = 1. From the known properties of the Laguerre 
polynomials one can show that in the cycle k there are k components of the Loschmidt echo 
(and the number of zeroes of the partial amplitude a^(t) is also k). The integral intensity 
of the echo within the cycle k is 2/T and independent of k, whereas the width of the echo 
(related to oscillation region for the Laguerre polynomials lOj) is ~ 4k. 

Because probability flow between the system and the reservoir states is determined by 
a s = —iCJ2n a n{t), one can describe the system functional space by the canonical variables 
a s (t) and a s (t). For a few initial cycle numbers, the trajectories are almost periodic in the 
weak coupling limit, and fill out densely the available phase space upon increasing of the 
coupling constant (see Fig. 2). For the higher cycles the trajectories become more and more 
complex and tangled, that is considered as a sign of chaotic behavior in classical dynamics 
[12I ]. 1^ | . The entanglement of the trajectories occurs when the echo components of the 
mixing for the overlapping cycles. To characterize this phenomenom one can introduce 
the critical cycle number fcW = ti 2 C 2 when the width of the oscillations for the Laguerre 
polynomials in the amplitudes a[ k \t) is equal to the cycle k period. For the k > kjp the 
cycles are overlapped. 

In the weak coupling limit upon increasing of the T, the Loschmidt echo intensity in- 
creases, while the coherent component of the amplitude decreases. We deduce a rough 
estimate of the critical cycle number k^') from the condition that the total duration of the 
cycles with the decay rate T = ttC 2 is of the order of the oscillation period 2tt/C. The 
criterion reads as 

k (2) „ 2C -3 _ 
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We conclude that there are three regimes of the time evolution which are shown in the 
Fig. 3. In the weak coupling limit (I) the coherent oscillations govern the dynamics, for the 
strong coupling region (II) - multicomponent Loschmidt echo takes place, and in the region 
(III) - the cycle overlap occurs. The crossover line (i) - (II) corresponds to the condition 
T = 1 also known as irreversibility criterion for the continuous reservoir spectrum (14J ]. In 
our case (discrete spectrum reservoir) when there are the recurrence cycles, this criterion 
means transition from coherent to incoherent behavior 15j |. In the region (I) the dynamics 
is ergodic, whereas in the region (II) it is not, since many rationally independent frequencies 
contribute to the time evolution. Note that the (I) - (II) crossover is similar to the known 

n 

in classical dynamics [13J ergodic - mixing transitions. 

The number of zeroes in the strong cycle overlapping condition (i.e., for k 2> fc£ ) is oc k, 
like for the non-overlapping cycles, however, only the small fraction (J k^ jk) of these zeroes 
belongs to the given cycle k. All other zeroes come from the previous cycles. Therefore the 
interval between the nearest zeroes in the cycle k is determined by a large number of the 
previous cycles. If in the cycle k the neighboring in time zeroes come from the previous cycles 
k', k", then in the next cycle k + 1 the interval between these zeroes increases. Moreover, 
in the cycle k + 1 at least one more new zero appears between those. These phenomena 
(mixing of zeroes and decreasing of the time interval between the neighboring zeroes) could 
lead to chaotic long time behavior at any finite accuracy of time or energy measurements (or 
by other words at any finite coarse graining of the system). Indeed in such conditions one 
may not restore uniquely the wave function from the time dependent amplitudes a s (t). To 
do it one has to know more and more precisely the widely oscillating function. At any finite 
accuracy one can always find the cycle number k c when dynamics of the system for k > k c 
can not be determined uniquely from the eigen state spectrum, and in this sense long time 
system evolution becomes chaotic. It is worth noting a similarity of this phenomenom with 
quantum mechanical uncertainty due to quantum object - measurement device interaction 



We do believe that our main results (multicomponent Loschmidt echo, three different 
regimes of the time evolution, cycle mixing) are valid at least qualitatively for a generic 
situation of any simple system coupled to a discrete spectrum reservoir, because in fact 
the results do not depend essentially on the details of the spectrum, but only on its average 
characteristics (the similar phenomenom is known for nuclear reactions [8|]). Note to the same 
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point that modern femtosecond spectroscopy methods (see e.g., |l| - |4j) indeed demonstrate 
(in a qualitative agreement with our consideration) remarkably different types of behaviors 
(exponential decay and complicated oscillation) of relatively close initially excited states. 
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Figure captions 

1. Fig. 1 

The system time evolution: (a) for C 2 = 0.3; (b) for C 2 = 1.0. 

2. Fig. 2 

The system trajectories for the overlapping recurrence cycles (C 2 = 1, [t/2ir] = 30). 

3. Fig. 3 

The phase diagram with three dynamical regimes: (I) - coherent oscillations; (II) - inco- 
herent evolution; (III) - mixing behavior. 
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